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Abstract 



The effect of an octupole term in a quadrupole deformed single particle poten- 
tial is studied from the classical and quantum mechanical view point. Whereas the 
^ ■ problem is nonintegrable, the quantum mechanical spectrum nevertheless shows some 

^ i shell structure in the superdeformed prolate case for particular, yet fairly large oc- 

tupole strengths; for spherical or oblate deformation the shell structure disappears. 
This result is associated with classical periodic orbits which are found by employing 
the removal of resonances method; this approximation method allows determination 
of the shape of the orbit and of the approximate octupole coupling strength for which 
it occurs. The validity of the method is confirmed by solving numerically the classical 
equations of motion. The quantum mechanical shell structure is analysed using the 
particle number dependence of the fluctuating part of the total energy. In accordance 
with the classical result this dependence turns out to be very similar for a superde- 
formed prolate potential plus octupole term and a hyperdeformed prolate potential 
without octupole term. In this way the shell structure is explained at least for some 
few hundred levels. The Fourier transform of the level density further corroborates 
these findings. 
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1 Introduction 



The shell structure is one of the most important quantum features of finite Fermi systems. 
The quantisation of a system of Fermions moving in a common potential leads to a bunching 
of levels in the single-particle spectrum, known as shells. When the levels of a bunch 
are filled the system is particularly stable; an additional particle fills a level in the next 
bunch at considerably higher energy and therefore produces a less stable system. In atoms 
and nuclei the shell effects are discerned by the deviations of the binding energies from 
the smooth variation obtained from the liquid-drop model or as oscillations in the radial 
density distribution as a function of particle number and of deformation [|^. As it was 
pointed out in the shells are a global phenomena not only for spherical systems. When 
a spherical shell is only partially filled, the higher degeneracy in the single-particle spectrum 
of a deformed mean field can lead to a breaking of spherical symmetry, and it therefore gives 
rise to a deformed equilibrium shape. 

A deep understanding of shell structure phenomena in terms of classical trajectories has 
been achieved by Balian and Bloch based on the periodic orbit theory of Gutzwiller ||^. 
According to the semiclassical theory the frequencies in the level density oscillations of 
single-particle spectra of nuclei are determined by the corresponding periods of classical 
closed orbits. The short periodic orbits give the major contribution to the gross shell 
structure ||, 0, 0- Depending on the particular mean field potential a deviation from 
spherical symmetry can lead to chaotic motion in the corresponding classical problem, and 
the shell structure of the corresponding quantum spectrum is affected or even destroyed 
depending on the degree of chaos |TD|, |TT|. 



Atomic clusters are another area which can serve as testing ground for these ideas. Recent 
experimental results on metallic clusters reporting abundance variations in mass spectra, 
ionisation potentials, static polarisabilities and collective giant dipole resonances, barrier 
shapes and fragmentation provide us with striking manifestations of shell structure effects 



related to a quantised motion of the valence electrons [|1^, |T3|. The correspondence of the 
electronic shell structure in spherical clusters to the closing of major quantal shells [0, 
caused considerable interest in using nuclear shell model type calculations for the description 
of metallic clusters fl^, It turns out that phenomenological potentials used 

traditionally in nuclear physics serve a purpose similar to those obtained within the Kohn- 



Sham density-functional method [21] if the relevant parameters are adjusted appropriately. 
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Typical potentials are the Woods-Saxon (and its various modifications) and the modified 
Nilsson potential without spin-orbit term. The considerable lowering of the computational 
time due to their simple analytical form renders an analysis of the stability of large metallic 
clusters feasible. Naturally, the shell numbers have to be larger than the ones used in nuclear 
physics in accordance with the larger number of valence electrons considered for metallic 
clusters. Recently, the semiclassical analysis has been successfully applied to explain shell 



structure produced by valence electrons p2| , |23| , Q in different metals. As it was mentioned 
before, obviously, shell structure in the quantum mechanical spectrum is associated with 
periodic orbits in the corresponding classical problem. Moreover, following work by Balian 
and Bloch Nishioka, Hansen and Mottelson [^] predicted the supershell, which is a 
beating pattern, when the regular oscillations refiecting the main shells are enveloped by 
an oscillating amplitude of lower frequency. As was argued in the case of a spherical cavity 
0, the beating pattern is caused by the interference between the two most important 
classical trajectories which are characterised by orbits of a triangular and square shape. For 
mesoscopic objects like clusters, deformations, i.e. deviations from spherical symmetry of the 
potential, are as important as in the nuclear physics context (for example, the smaller peaks 
in the cluster abundances correspond to unfilled shells or deformed shapes of clusters p5|). 
Direct information about the deformation of clusters can be obtained from an analysis of a 
splitting of the plasmon resonance or from a more detailed investigation of mass abundance. 

Based on a simple analysis of the geometry of closed classical orbits, Mottelson |2^, |1| 
predicted a tendency for a large system (heavy nuclei) with a 2:1 quadrupole deformation 
to exploit the octupole degree of freedom. For nuclei this idea was supported by Nils- 
son+Strutinsky shell-structure energy calculations [^], by a systematic study of octupole 
correlations within a microscopic approach |2D[ and by a qualitative analysis of the 
dynamical symmetry of the anisotropic harmonic oscillator [3D|. Hamamoto et al |HT 



m- 



vestigated the importance of nonaxial octupole deformation for many-body system within 
a harmonic oscillator model. More recently, authors in [jl8l used a deformed Woods-Saxon 
potential, including quadrupole, octupole and hexadecupole terms, to calculate the ground- 
state deformation of sodium clusters. 

Inclusion of an octupole term in addition to a quadrupole term renders the classical 
single particle motion nonintegrable. In fact, the system turns out to be chaotic. In a recent 
investigation the study of classical motion in a cavity with oscillating walls of even and 
odd higher order multipoles has led to interesting conclusions about elastic versus dissipative 
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behaviour of a noninteracting gas depending on the integrabihty or nonintegrabihty of the 
equations of motion. 

In the present paper we thoroughly investigate a simphfied model. We leave out terms, 
which, albeit physically important, are prone to blur the analysis when the interest is fo- 
cussed on the essentials, in particular a distinction between orderly and chaotic motion in 
many body systems, like nuclei and metallic clusters. Since our interest is directed not only 
towards the classical but also the corresponding quantum mechanical motion, we leave out 
the spin-orbit term and the (Z)^-term present in the Nilsson model so as to render as closely 
as possible the analogy between the classical and quantum cases. An analysis of classi- 
cal trajectories in a quadrupole deformed harmonic oscillator has shown that the spin orbit 
term leads to chaotic motion Likewise, the (/)^-term gives rise to chaotic behaviour |^ 



even without the octupole term. Despite and because of the simplifications, such a model 
allows an understanding of the main features of shell structure effects in super (hyper) de- 



formed nuclei [0, and in deformed clusters (see, for example, |jT2|, ^ and references 
quoted therein). Recent experimental data of superdeformed K-isomers in nuclei and 



electronic shell structure effects in metallic clusters [|T2[ seem to indicate the manifestation 
of oblate deformation in nuclei and small clusters. In the present paper we investigate the 
effect of the octupole term for prolate and oblate deformation; the latter case has not been 
dealt with in 0. 

One major result of the present paper is a classical analysis which demonstrates that 
the prolate case including octupole deformation is still quasi-integrable . In this way the 



quantum mechanical shell structure found earlier [TTI is given a proper theoretical foun- 
dation. One may speculate that this is the reason why for large quadrupole deformation 
prolate-octupole nuclei (clusters) are more stable than the oblate-octupole ones. For metallic 
clusters finite temperature should be considered 0; however, our interest is focussed on 
shell structure whose character is unaffected by temperature except for the amplitude. The 
detailed analysis of magic numbers, their changes under variation of the octupole strength 
and the possibility of supershell structure induced by the octupole term is a further major 
contribution of this paper. 
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2 The model 



We investigate the classical and quantum mechanical motion in the potential 

v(...)^f.v4 + A?£l;||!,. (1) 

For 6 > 1 (6 < 1) this is a quadrupole deformed harmonic oscillator of prolate (oblate) shape 
with an additional octupole term; in fact the term proportional to A is r^P3(cos^) with P3 
the third order Legendre polynomial. We use cylindrical coordinates z and q = y/x'^ + y'^. 
For A 7^ this is a two degrees of freedom system which is non-integrable. 

The exact classical orbits are obtained numerically by integrating the equations of mo- 
tion, viz. 

z = -lvig,z). (2) 

It is convenient to use the energy as one initial condition; because the potential obeys the 
scaling law V{'yg,'-yz) = 'y'^V{g,z) one value of the energy yields orbits of all energies. The 
equation for g would contain the term p'^/ g^ on the right hand side of Eq.(^ if a finite value 
of the (conserved) angular momentum p,^ would be considered. However, as discussed in 
the following section, = is sufficient for all classical aspects relevant in this paper. 

Our choice of the octupole term ensures that we have a genuine bound state problem 
for A < Acrit- Here Acrit is defined to be the value for which the potential no longer binds, 
and for |A| > Acrit the potential tends to —00 along one or two directions in the p-z-plane. 
The direction and the value of Acrit depend on the quadrupole deformation b. For prolate 
nuclei (6 > 1), Acrit = l/(2&^) and the potential opens its valley along the positive (negative) 
2;-direction for negative (positive) A. For oblate nuclei (6 < 1) the other possible direction 
along the line g = sign{X)j3z with (3 ~ 0.4, is of increasing importance. At the value b ^ 0.58 
valleys along the two directions ^ = and g = OAz open simultaneously for Acrit ~ 1.5 while 
for a still smaller value of b, say for b = 0.5, the valley along the direction g = OAz opens 
for Acrit ~ 1-64 while the one along the 2;-axis now opens for a larger value of A. Analytic 
expressions for Acrit and the direction j3 as functions of b are given in Appendix A. 
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3 Classical perturbative treatment 



To understand the occurrence of shell structure in the corresponding quantum mechanical 
problem we present a classical analysis in which it becomes obvious that, for the prolate 
case, the problem is approximately equivalent to an integrable problem even for finite values 
of A. In fact it is shown in what follows that, within the approximation used, the motion 
in the z and ^^-direction becomes uncoupled whereby the motion in the z-direction depends 
on A. In this way frequencies in either direction can be clearly defined with the effect that 
the corresponding winding number, i.e. the ratio uOg/ujz-, becomes a simple function of A. 
This provides the basic mechanism for the occurrence of shell structures for the quantum 
mechanical problem; in fact it provides classical evidence for the quantum mechanical finding 
that the spectrum for specific values of A is nearly equivalent to a spectrum with A = but 
with a larger value of the quadrupole parameter h. We stress that the agreement between 
the perturbative results and the exact orbits is excellent. 

The method employed is based on the 'removal of resonances' approach which amounts 
to a particular averaging procedure used in secular perturbation theory |^ . The complete 
Hamilton function written in terms of the angle and action variables of the unperturbed 
problem (see Appendix B) reads 

/T , It , . y5J2sina^(26J2sin^a2 - SJgSin^ag) 
H[Jg,J,,ag,a,)=uj[Jg + -J, + \ ). (3) 

" J bJz sin az + JoSin an 



The frequencies Uz = §j- and Ug = are given by 

^ uj 6sin 02(46^ J| sin^ + GfoJ^Jr, sin^ a„ sin^ — 3J^ sin 4a, 



b 2\/bJz{bJz sin^ + Jg sin^ a^) 2 

- sinazSin^ ag^/Wz{%bJzSin^ az + ^JgSiv? ag). , 

Ug{J,a) = u{l-X , J . 2 O )• 

2[bJzSm a^ + J^sm ag)2 

It is clear from these expressions that, for small values of A, the ratio ojg/uoz is essentially 
determined by b. For prolate superdeformation, when 6 = 2 or larger, one can expect 
that averaging over the fast moving angle Ug of the Hamilton function should yield a good 
approximation. Since the Hamilton function is periodic in both angles, this approximation 
amounts to keeping the zero order term of its Fourier expansion in the fast moving angle 
[|37| . Since integration of the Hamilton function over Ug can be done analytically, this 
approximation appears particularly attractive. 
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After averaging and then rewriting the action variables and the remaining angle in 
terms of the original momentum and coordinate values we obtain from Eq.@ 

where = 2J^^^/{muj) = + p^^/ {mujY which is a constant within the approximation. 
Here F21 and K are the hypergeometric function and the first elliptic integral, respectively. 

With Eq.@ we have effectively reduced the full problem to an integrable problem as the 
motion in the two coordinates is uncoupled. The motion in the 2;-coordinate is determined by 
an effective potential which now depends also on the action Jg in addition to its dependence 
on A. Note that the effective potential originating from the octupole term depends only 
on z/^, ^ ~ \J~Tq] this means that, depending on the action residing in the ^)-motion, the 
amplitudes of the vibrations in the 2;-motion are increased or decreased. The potential of 
the 2-motion is no longer harmonic. Because of its dependence on A and Jg we now expect 
ratios ujg/uOz which are different from the unperturbed value h and which depend on A and 
Jg. Also note that the large z behaviour coincides with that of the complete problem which 
means that Acnt = l/(2&^) still holds. In Fig.(l) we illustrate the effective potential. 

A few comments are in order with regard to non-zero angular momentum. The additional 
term p^/ in the Hamilton function does not alter our procedure nor the results based 
upon it. In fact, the f?- motion remains unperturbed including the centrifugal term, and 
for the harmonic oscillator the frequency ujg remains therefore equal to its zero angular 
momentum value. This is significant in the quantum mechanical context as the ratios ujg/ujz 
are independent of the angular momentum just like in an unperturbed harmonic oscillator. 
It is by this mechanism that we can explain pronounced shell structure for A 7^ 0, i.e. shells 
containing all values of angular momentum which are allowed quantum mechanically. We 
return to this point in the subsequent section. 

The approximation procedure seems to be useful only for the super- and hyperdeformed 
prolate case. In the superdeformed oblate case (b=l/2) one would, at first glance, consider 
averaging over a^. However, the zeroth order term vanishes now as the octupole term is an 
odd function of smaz- The first order term in the Fourier series no longer decouples the 
coordinates q and z so that no gain is made. For the spherical case (6 = 1) a canonical 
transformation could in principle enforce one of the (transformed) frequencies to be faster 
than the other. However, this procedure again is not useful since the slow frequency turns 
out to be a complicated function of time, in other words, the ratios ojg/ujz fail to be con- 
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stant rational numbers. In simple terms, this discussion shows that only the fairly strongly 
deformed prolate case remains close to integrability when the octupole term is switched 
on. This is the classical explanation for the existence of shell structure in the quantum 



spectrum which occurs only for the prolate case |]rT|. As we see below the spherical and 
the oblate case soon become chaotic classically when the octupole is switched on, which 
prevents occurrence of shell structure in the corresponding quantum spectrum. 

The frequency ujz as determined by the approximation implied by Eq.(^ is given by 
ujz = 27t/T where 

T = V2^ / , (6) 

^E, - U{z) 



(7) 



with Ez = E — Eg = E — miJ^^ jl and Ui^z) is taken from Eq.(^ and reads 



U{z) = 



The period T can be evaluated numerically. In Fig. (2) we display the ratio uJq/uJz as a 
function of A for a few parameter values ^. Starting at its unperturbed value 6 = 2 for 
A = 0, it turns out that the ratio of the frequencies is a monotonously increasing function 
of A for < A < Acrit for a large range of ^-values, for which it is virtually independent 
of ^; only when approaches its upper limit 2E/{muj'^) can a notable ^-dependence be 
discerned. We note that the integral in Eq.(|D can be evaluated analytically for ^ = 0; in 
this case the ratio reads 

' + ' ) (8) 



2 ^1 + A/A„i, \/l-A/A„,.' 



which is the top curve drawn in Fig. (2). Simple, i.e. short periodic orbits occur when the 
frequency ratio becomes a rational number p/q with small values of p and q. From Fig. (2) 
we read off that around A = 0.66Acrit the ratio is 5:2 and for A = O.SAcrit it is 3:1. In Figs. (3) 
we display the most important periodic orbits. If the approximation is meaningful this 
finding implies that for these particular values of A all orbits (except the ones with Ez > 0) 
should be periodic with the respective winding numbers. The consequence for the quantum 
spectrum is clear: the parameters 6 = 2 and A = 0.66Acrit should yield a spectrum similar to 
that for the parameters 6 = 5/2 and A = 0; likewise, 6 = 2 and A = O.SAcrit should be nearly 
equivalent to 6 = 3 and A = 0. In the following section we see a convincing confirmation of 
this classical prediction of quantum shell structure. 
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One expects a similar pattern when starting from the outset, say, with b = 5/2, since 
the onset of classical chaos seems to be delayed the larger the value of b. In fact, one finds 
appropriate A-values for which the winding numbers become 3:1 (for A = 0.63Acrit) and 7:2 
(for A = 0.76). Starting with b = 3 one finds 7:2 and even 4:1, i.e. we have the the "chain 
of shells" 2:1-^5:2-^3:1^7:2-^4:1. ..We mention here the 7:2 orbit to 
demonstrate the principle, but we have left out the 7:3 in the previous paragraph although 
it is there between the 2:1 and 5:2. As a matter of fact we have left out all ratios p/q with 
p > 7 because they do not feature as prominently in the quantum spectrum as the shortest 
orbits with p < 7. But we can draw the conclusion that an octupole term in a prolate 
deformed oscillator produces essentially a pure but more deformed prolate oscillator. This 
explains the occurrence of shell structure in the corresponding quantum spectrum at the 
particular A-values for which the effective value of b assumes a rational number with small 
value of the numerator and denominator. 

We have tested the validity of the approximation procedure by looking at a large variety 
of exact orbits obtained from solving Eq.(^ numerically. In Fig. (4a) surfaces of section 
are displayed for A = 0.63Acrit for a number of initial conditions. The slight discrepancy 
between the values 0.66 found from Fig. (2) and 0.63 used in Fig. (4a) reflects the difference 
between the approximation and the exact situation. Like in Fig. (2) the actual ratio when 
determined numerically is an increasing function of A. It turns out that the actual ratio is 
slightly larger than 5:2 for A = 0.66Acrit and is closest to 5:2 for A = 0.63Acrit- What meets 
the eye when looking at Fig. (4a) are the large stability islands around the 5:2 orbit. Note 
that at the periphery of the surface which is chosen to be the g = plane we have = 
while attains its maximum value towards the centre. In this way we sample the whole 
range. Of course we do not obtain an exact 5:2 orbit for a large range of .^-values as the 
approximation would suggest; we do get the one in the centre of the islands. However, while 
most orbits are quasiperiodic and produce tori, their winding numbers are all very close to 
5:2, for instance 70:27 or 70:29. This applies to most of phase space, deviations are noticeable 
only towards the centre of the section. In accordance with Fig. (2) there we get orbits with 
smaller winding numbers such as 7:3 and even 3:2. However, more than 85% of phase space is 
dominated by orbits with winding numbers very close to 5:2. The analogy of the parameter 
set 6 = 2, A = 0.63Acrit with b = 5/2, A ~ goes even further in that not only the shape of 
the trajectory of the actual 5:2 orbit is the same but also that of the unstable orbit associated 
with the intersection points of the separatrix (a separatrix similar to the one displayed in 
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Fig. (4a) is produced by the parameter set b = 5/2, A = 0.05Acrit)- Mutatis mutandis similar 
statements hold for the comparison of the parameter set 6 = 2, A/Acrit = 0.76 (see Fig. (4b)) 
and 6 = 3, A ~ 0. Note the small region where the onset of chaotic motion can be discerned; 
by and large the motion still appears close to integrabilty although A is fairly close to Acrit- 
For comparison we display surfaces of sections in Fig. (5) for the spherical and oblate case. 
Even though the diagrams represent results for smaller values of A/ Acrit the overall onset of 
chaos is obvious. Here we note that for the special value 6 ^ 0.58 discussed in section 2 the 
onset of chaos occurs for an even smaller value of the octupole strength, i.e. A/Acrit ^ 0.25. 

4 Quantum mechanical results 

The quantum mechanical treatment is straight-forward in principle. In the spirit of previous 
work we use for the full problem which is of the form Hq + XHi a representation with 
Hq diagonal. The basis chosen is referred ^ to as the basis using the asymptotic quantum 
numbers n±, and A where n± = n+ + and A = — n^. Here the quantum numbers 
n+ and n_ are the eigenvalues of (^4.)"''^+ and {A^YA^, respectively, where, in terms of the 
usual boson operators ax and ay, we use A± = (a^ =F io.y) / V^- For a fixed value of A (which 
is the quantised analogue of the classical p^) this leaves two quantum numbers (reflecting 
the two degrees of freedom) to enumerate the rows and columns of the matrix problem. For 
A = the diagonal entries of Hq are thus £nj_,n:, — + 1 + (n^ + l/2)/6). The matrix 

elements of Hi are obtained from those of z ~ (a J + a^) and g"^ ~ {A-^-{A^y + A^{A^y + 
A+A_ + {A^y (A^Y). To get the matrix elements of + z'^ in a numerically consistent 

way, we first calculate the matrix elements Sm,n of + from which the inverse square 
root is obtained using 

5-1/2 ^ jj . jj-i/2 . jj] ^Yiem D = ■ S ■ U is the diagonal form 
of the positive definite matrix S and U is the orthogonal matrix which diagonalises S. To 
ensure also numerically that S has only positive eigenvalues it is important that the matrix 
for z"^ is obtained by squaring z and not by evaluating analytically the matrix elements from 
{al + a^y-, inconsistencies are otherwise introduced due to truncation. In this way, we also 
ensure that the truncated matrices S~^^'^ and the representation of 2z^ — 3zg^ commute and 
that their product is a symmetric matrix. The effect of truncation was tested by looking 
at the variation of the lower end of the spectrum when the dimension of the matrices was 
increased. There is certainly a dependence on 6 and A. For A < 0.9Acrit and 0.5 < 6 < 2 
the variation was less than 1% for the first 300 levels obtained from 1600x1600 dimensional 
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matrices. 

In Fig. (6a) we illustrate the spectrum so obtained as a function of A for b = 2. The 
shell structure at about A = 0.63Acrit and A = 0.76Acrit is clearly discernible. While Fig. (6a) 
presents the spectrum for A = we illustrate in Fig. (6b) the whole spectrum which is 
a superposition of all possible A-values. Here the shell structure is even more pronounced 
which is expected since, according to the discussion in section 3, the orbits will have the same 
winding numbers independent of the angular momentum; as a consequence, the quantum 
spectra will have shell structures similar to an unperturbed oscillator. This is a prototype 
case for employing Strutinsky's procedure Q to exhibit the magic numbers associated 
with the shells. However, while we follow the basic idea of Strutinsky's method, we use 
a somewhat different approach which appears more suitable in our case; this is presently 
discussed. 

It is clear that for the plain quadrupole deformed oscillator, supershell structure in the 



sense of p2[ cannot occur since there is basically only one orbit for each coordinate; 
it is the sinusoidal motion characterised by the respective frequency associated with that 
coordinate. The consequences for the quantum spectrum are well known [||]: one obtains the 
Fibonacci numbers as degeneracies in the equidistant spectrum. If the frequencies ratios 



are simple rationals, there are systematic repetitions of the same degeneracies |^0|, ^ 
This is in contrast to, say, the motion in a cavity or in a Woods-Saxon potential where 
basically different classical orbits such as the triangle, the square, the pentagon and so forth 
exist thus giving rise to interference of different shells in the quantum spectrum. One of 
the questions addressed in this paper is the possible occurrence of supershell effects when 
a quadrupole deformed harmonic oscillator is perturbed by an anharmonic term, i.e. the 
octupole deformation. 

To analyse the quantum spectrum we proceed in the orthodox fashion in that the total 
energy EtotiN,\) = Y^i^ is approximated by a smooth average function S{z,X) and 
the fluctuating difference 

SE{N,X) = E,,,{N,X)-S{N,X) (9) 

is then subjected to further investigation. The finding of a suitable form for the average 
function S{z, A) is facilitated in our case as it is well known that the leading term of 
EtotiN, X), as a function of N, is proportional to N'^^^. We determine the five constants 
ao(A), . . . , a4(A) in S{z,X) = J2k=o'^kiX)z''^^ by a least square fit which turns out to be 
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perfectly satisfactory for all values of < A < Acrit- 

In Fig. (7) the fluctuating part SE{N,X) is presented as a contour plot. The diagram 
refers to 6 = 2 and displays the ranges 100 < < 700 and < A < Acrit- For A = (the 
bottom horizontal line) we clearly discern the shell structure of the plain deformed oscillator. 
Note that the sharp minima (dark shadowing) occur at the positions where a shell is 
closed; the distances arc proportional to N^. When A is switched on the shell structure 
persists to a great extent; only when A approaches its critical value (top horizontal line) the 
structure begins to be washed out. There are local minima discernible at A/Acrit ~ 0.76. 
This is a reflection of the enhanced shell structure discussed in the previous section. Below 
we pursue this point in more detail. The finer structure and in particular the detailed shell 
structure is obtained from the second derivative of 5E{N, A) to which we turn next. 

In Figs.(8) we have plotted g{E) = 6E{N + 1, X) + 6E{N -1, X) -26E{N, X) as a function 
of N^^^ for a few characteristic values of A. The peaks of the plots represent the magic num- 
bers which characterise the shells, and the heights of the peaks reflect the energy distance 
from one shell to the next. In Fig.(8a) the essentially unperturbed result (A/Acrit = 0.15) 
is presented for demonstration. The second row displays the results for the particular val- 
ues of A for which the winding numbers are 5:2 (A/Acrit = 0.63) and 3:1 (A/Acrit = 0.76), 
respectively. The magic numbers and the heights of the peaks agree well with those which 
are obtained from the unperturbed (A = 0) quadrupole deformed oscillators (third row in 
Fig.(8)) with b = 5/2 and 6 = 3, respectively, at least for N < 700, A^^s < g.gg. xhe 
agreement extends in particular to the respective occupation numbers, i.e. the degeneracies; 
of course, the heights of the peaks do not show the same regularity as the correspond- 
ing unperturbed problem; nevertheless, even for the heights an overall agreement prevails 
when comparing with Figs. (Be) and (8f) where the respective unperturbed structures are 
displayed. For higher values of N we do get deviations which reflect upon the fact that 
the system is nonintegrable and cannot give complete order in all its results. While the 
agreement for lower values of N was to be expected from the discussion in the previous 
section, the extent of the agreement is rather remarkable, especially for A ~ 0.76Acrit where 
an astoundingly clean shell structure reoccurs after it partially disappeared for a somewhat 
smaller value of A. It is this recurrence of shell structure which gives rise to the local minima 
in Fig. (7) as pointed out above. 

Of particular interest is Fig. (8b) which refers to the intermediate value A = 0.7Acrit 
where the genuinely different orbits with winding numbers 5:2 and 3:1 coexist. The long 
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wave length fluctuation could well be interpreted as a supershell structure, even though 
that it is not as clearly pronounced as in a more transparent integrable case Yet the 
difference to Fig.(8d) which refers to a larger value of A is striking. 

In Figs. (9) we display the square of the modulus of the Fourier transform of the level 
density, i.e. the function 

= |^e'^"f . (10) 

n 

The spectrum is taken at A = 0.63Acrit (Fig-(9a)) and at A = 0.76Acrit (Fig(9b)), both 
spectra refer to A = only. The pronounced peaks can be directly associated with the 
periods of the classical 5:2 and 3:1 orbit, respectively, the periods obtained from Figs. (9) 
are in perfect agreement with those of the corresponding classical orbits which are found 
numerically by integrating Eqs.(0). This is a beautiful demonstration of Gutzwiller's trace 
formula . As expected the frequencies deviate considerably from the unperturbed values, 
i.e. from the frequency associated with 6 = 2, A = 0, but also from the frequencies associated 
with h = 5/2 or 6 = 3. In units of the unperturbed value (6 = 2) we find T5.2 = 1.2 and 
T3.1 = 1.4; the values are larger than unity in accordance with Fig. (2). 

Again we stress that the high degree of order which prevails in the superdeformed prolate 
case when the octupole term is turned on, does not exist in the corresponding oblate, in 
fact, not even in the spherical case. There, chaotic behaviour becomes manifest for much 
lesser octupole strength, which results in a complete disappearance of shell structure in the 
quantum spectrum. 



5 Summary 

The main motivation of our work is the understanding of the nature of shell effects of many- 
body systems like nuclei or metallic cluster at large quadrupole/octupole deformation. It 
is the shell structure that is responsible for the existence of superdeformed nuclei which 
nowadays is a very broad research subject. Metallic clusters provide us with an additional 
challenge; this topic is attracting more and more researchers who have traditionally been 
working in nuclear structure. In spite of the different character of the interaction in nuclei 
(predominantly attractive forces between nucleons) and metallic clusters (repulsive forces 
between valence electrons), the remarkable similarity in the manifestation of shell structure 
in experimental data allows analysis of both systems within the same quantum mechani- 
cal model. We have investigated the shell structure produced by a quadrupole deformed 
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harmonic oscillator with an octupole term. Contrary to the case of spherical potentials 
our Hamiltonian is nonintegrable. Using the method based on the 'removal of resonances' 
approach we established a connection between shell structures in the quantum me- 
chanical spectrum and periodic orbits in the corresponding classical problem. In the prolate 
case the classical problem can effectively be reduced to an integrable set of equations, since 
the motion in the two coordinates becomes uncoupled. This is in contrast to the spherical 
and oblate cases, where the motion becomes chaotic when the octupole term is switched 
on. In the prolate case, at particular values of the strength parameter A the octupole term 
produces a motion which resembles to a great extent that of a pure but more enhanced 
prolate oscillator. This provides the classical explanation for the existence of quantum shell 
structure for prolate/octupole deformed system within the model |TT|. At small values of 
A the most important short orbit is the Lissajou figure with the shape of an 'eight' with 
winding number 2:1 (top of Fig. (3)). 

The fluctuating part of the energy has been extracted from the quantum mechanical 
single particle spectrum using our procedure which is a variant of the method discussed 
in 0. We have found remarkable agreement between manifestations of shell structure for 
the same values of A/Acrit (Fig- 7) and the ones which give rise to stability islands in the 
Poincare surfaces of sections relating to the classical short orbits with winding numbers 2:1 
(A ^ 0), 5:2 (A/Acrit = 0.63), 3:1 (A/Acnt = 0.76) (Fig.(4)). Note that the orbit closest to 



the triangle orbit, which is one of the most important orbits in a cavity occurs for 

the quadrupole/octupole deformed harmonic oscillator only at a large value of the octupole 
strength parameter A; the orbit has the winding number 3:2 and is associated with the 
three islands in the centre of Fig. (4b), but it does not feature in the quantum mechanical 
spectrum. In the intermediate case (A/Acrit = 0.7) when the orbits with winding numbers 
5:2 and 3:1 coexist the long wave length fluctuation could be interpreted as a supershell 
structure. But due to the narrowness of these orbits the supershell structure is not as 
clearly pronounced as in a more transparent integrable case |2^ . 

In this paper we placed our emphasis on the superdeformed case which is the most in- 
teresting situation when the octupole term is switched on. Shell structure is destroyed for 
smaller values of the octupole strength for lesser quadrupole deformation or even more so 
for the spherical case. Our analysis implies that the shell structure favours the superde- 
formed prolate/octupole in contrast to the oblate/octupole case. The model is only a crude 
approximation of a realistic situation, and a mean field without consideration of correlations 
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(interactions between particles) may not justify a statement as strong as this. Only future 
experimental data from the newly built facility like EUROGAM and GAMMASPHERE can 
test our conclusion with regard to nuclei. For clusters one could expect that the process of 
cluster aggregation and cluster evaporation will follow valleys of the potential energy surface, 
which is the combined effect of electronic shell structure and ionic lattice. The electronic 
shell effects should be rather important at relatively small temperature and not for too large 
number of particles. In the framework of our model the electronic shell structure becomes 
less pronounced for > 700, since higher up in the spectrum the chaotic behaviour may 
well interfere with the search for shell structure |]3^ . Alternatively, if shell structure can be 
observed for larger numbers, our conclusion would be that weaker deformation must prevail. 
More experimental evidence is needed to assess the situation; such evidence is expected also 
to shed light on the question of whether clusters provide further insight on the fascinating 
question of quantum behaviour of classically chaotic systems. 

Finally we comment on the (/)^-term in the Nilsson model, which plays a positive role as 



a phenomenological term in nuclear physics. From our results one should infer that this 
term destroys shell structure, since it produces chaos in the classical analogous (quadrupole 
deformed) case, while the experimental data indicate shells for clusters. This is different 
from nuclear physics where chaos does not show, since one is not interested in too many 
levels. Further investigations have to provide clarity on this aspect. Work to this effect is 
in progress. 



RGN gratefully acknowledges financial support from DGICYT of Spain. 



16 



A Determination of /3 and Acrit 



To study the behaviour of the potential along the line z — Pgwe consider the expression 

2l/te.)/(,W) = l + ^ + A^^. 

The zero of the derivative gives the maximum/minimum of the leading behaviour of the 
potential along the line z = (3q. From this we find the relation 



Solving for Acrit we obtain 



'^crit 



&2(4/34 + 6/32-3)' 

Inserting this for A into the previous equation yields the slope /3 as a function of h when 
A = Acrit, i-e. 

J_ I 9 v/3(3 + 26^)(l + 146^) ' 

^ + ^+ 2 + 62 

B Action and Angle Variables 

The relations between the coordinates and momenta and the actions Jj and angles ctj are 
given by 



2J, . 
Q — \ — - sm a. 



1 2b J, 



sm az 



Pg — J2mujjQC0sag 



2mujJg 
Pz = \ — cosa^. 
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Figure captions 



Fig.l The potential U{z) in units of ma;^/2 for a few values of A and an intermediate 
value of ^. Variation of ^ produces basically similar shapes. 

Fig.2 The ratio ioJiOz as a function of A for various values of ^. A genuine ^-dependence 
is discernible only for ^ > O.S^max- 

Fig.3 The major orbits with winding numbers 2:1 (A > 0), 5:2 (A/Acnt = 0.63) and 3:1 
(A/AeHt = 0.76). 

Fig.4 Surfaces of sections (6 = 2) of four orbits for A = 0.63Acrit (top) and of five orbits 
for A = 0.76Acrit (bottom). 

Fig.5 Sections of three orbits for 6 = 1 (top) and for h — 1/2 (bottom) for A = O.SAcrit- 

Fig.6 Spectrum (6 = 2) for A = (top) and all A (bottom) as a function of A. The 
energy unit is huj. Note that there are only avoided level crossings in the top figure which, 
however, are not resolved. 

Fig.7 Contours of 5E{N,X) with 100 < < 700 as abscissa and < A/Acrit < 1 
as ordinate. Dark areas represent minima. Volume conservation is taken into account in 
Figs. (7) and (8). 

Fig.8 Shell structure for fixed values of A and b — 2. The second derivative of 5E{N, A) 
is illustrated versus N^/^ for A/Acrit = 0.15,0.70,0.63 and 0.76 in (a), (b), (c) and (d), 
respectively. In (e) and (f) the unperturbed (A = 0) structure is displayed for h — 2.5 and 
6 = 3, respectively. 

Fig.9 Fourier transforms of the modulus square of the level density for A = 0, 6 = 2 for 
A/Acrit = 0.63 (top) and A/Acrit = 0.76 (bottom). 
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